Is Thick Brane Model Consistent with Recent Observations? 
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There exist many observational evidences implying the expansion of our universe is undergoing a 
late-time acceleration, the mechanism of this acceleration is yet unknown. In the so-called thick 
brane model this phenomena is attributed to the thickness of the brane along the extra dimension. 
In this study we mainly rely to the consistency of this model with most recent observational data re- 
lated to the background evolution. The new Supernova Type la (SNIa) Gold sample and Supernova 
Legacy Survey (SNLS) data, the position of the acoustic peak at the last scattering surface from the 
Wilkinson Microwave Anisotropy Prob (WMAP) observations and the baryon acoustic oscillation 
peak found in the Sloan Digital Sky Survey (SDSS) are used to constrain the free parameter of the 
thick codimension 1 brane model. To infer its consistency with age of our universe, we compare the 
age of old cosmological objects with what computed using the best fit values for the model param- 
eters. When the universe is matter dominated, w = 0, at 68% level of confidence, the combination 
of Gold sample SNIa, Cosmic Microwave Background (CMB) shift parameter and SDSS databases 



provides Q. 



0.3li°;<g, fic 



0.05 



+0.01 



-1.401q'2q, hence a spatially open universe with 



0.211q'q8- The same combination with SNLS supernova observation gives fl r , 



0.281 



flc = 0.037 



+0.003 



w r = -2.05 



+0.15 



consequently provides a spatially open universe = 0.1L 



0.02 • 
0.10 



These results obviously seem to be in contradiction with the most recent WMAP results indicating 
a flat universe. 

PACS numbers: 98.80.Es, 95.36.+X 



I. INTRODUCTION 

Recent observations of type la supernovas (SNIa) sug- 
gest that the expansion of the universe is accelerating 
[1-4]. As is well known all usual types of matter with 
positive pressure generate attractive forces, which decel- 
erate the expansion of the universe. Given this, a dark 
energy component with negative pressure was suggested 
to account for the invisible fuel that drives the current 
acceleration of the universe. There are a huge number of 
candidates for the dark energy component in the litera- 
ture (sec, e.g., [5-10] for recent reviews), such as a cosmo- 
logical constant A [6,11-13], an evolving scalar field (re- 
ferred to by some as quintessence: [14-26], the phantom 
energy, in which the sum of the pressure and energy den- 
sity is negative [19,27,28], the quintom model [17,29], the 
holographic dark energy [30], the Chaplygin gas [31,32], 
and the Cardassion model [29,33,34]. Another approach 
dealing with this problem is using the modified gravity 
by changing the Einstcin-Hilbcrt action. Some of models 
as 1/R and logarithmic models provide an acceleration 
for the universe at the present time [33,35,36]. 
In addition to the phcnomenological modifications of the 
action, the brane cosmology also implies modification 
on the dynamics of the universe. Many different Brane 
models have been presented in the recent years. In one 
of the best studied models, a codimension 1 thin brane 
with an infinite extra dimension has been investigated. 
This brane is embedded in a bulk where its curvature 
to be negative and its volume is finite [37]. In this sce- 



nario, gravitation is localized on a brane reproducing ef- 
fectively four-dimensional gravity at large distances due 
to the warp geometry of the spacctimc. Usually the 
brane is modeled as a distributional source in the energy- 
momentum tensor (EMT) of zero thickness, and in this 
case the cosmology has been obtained and analyzed in de- 
tail [38,39]. Recognizing the difficulty of handling thick 
walls within relativity, already early authors considered 
the idealization of a singular hypersurface as a thin wall 
and tried to formulate its dynamics within general rela- 
tivity [40] . The new era of intense interests in thin shells 
and walls began with the development of ideas related to 
phase transitions in early universe and the formation of 
topological defects. Again, mainly because of technical 
difficulties, strings and domain walls were assumed to be 
infinitcsimally thin [41]. 

Thereafter, interest in thin walls, or hypcrsurfaccs of dis- 
continuity, received an impetus from the cosmology of 
early universe. The formulation of dynamics of such sin- 
gular hypersurfaces was summed up in the modern ter- 
minology by Israel [42]. Within the Sen-Lanczos-Israel 
(SLI) formalism, thin shells are regarded as idealized zero 
thickness objects, with a (^-function singularity in their 
energy-momentum and Einstein tensors. 
In contrast to thin walls, thickness brings in new sub- 
tleties, depending on how the thickness is defined and 
handled. Early attempts to formulate thickness, being 
mainly motivated by the outcome of the idea of late 
phase transition in cosmology [43] , were concentrated on 
domain walls. Bonjour et al., studied a thick gravitat- 
ing domain wall for a general potential [44,45]. Using 



general analytical arguments they have shown that all 
nontrivial solutions fall into two categories: those inter- 
pretable as an isolated domain wall with a cosmological 
event horizon, and those which are pure false vacuum 
de Sitter solutions. Also they have analyzed the domain 
wall in the presence of a cosmological constant finding 
the two kinds of solutions, wall and de Sitter, even in 
the presence of a negative cosmological constant. Silve- 
ria [46] studied the dynamics of a spherical thick domain 
wall by appropriately defining an average radius for the 
wall. Widrow [47] used the Einstein-scalar equations for 
a static thick domain wall with planar symmetry. He 
then took the zero-thickness limit of his solution and 
showed that the orthogonal components of the energy- 
momentum tensor would vanish in that limit. Garfinklc 
and Gregory [48] presented a modification of the Israel 
thin shell equations by using an expansion of the cou- 
pled Einstein-scalar field equations describing the thick 
gravitating wall in powers of the thickness of the domain 
wall around the well-known solution of the hyperbolic 
tangent kink for a A(f> 4 wall and concluded that the ef- 
fect of thickness at first approximation was effectively to 
reduce the energy density of the wall compared to the 
thin case, leading to a faster collapse of a spherical wall 
in vacuum. Others [49] applied the expansion in the wall 
action and integrate it out perpendicular to the wall to 
show that the effective action for a thick domain wall in 
vacuum apart, from the usual Nambu term, consists of a 
contribution proportional to the induced Ricci curvature 
scalar. 

Study of thick branes in the string inspired context of 
cosmology began almost simultaneously with the study 
of thin branes, using different approaches. Although 
in brane cosmology the interest is in local behavior of 
gravity and the brane, most of the authors take a pla- 
nar brane for granted [50]. However, irrespective of the 
spacetime dimension and the motivation of having a wall 
or brane, as far as the geometry of the problem is con- 
cerned, most of the papers are based on a regular solu- 
tion of Einstein equations on a manifold with specified 
asymptotic behavior representing a localized scalar field 
[51]. Some authors use a smoothing or smearing mech- 
anism to modify the Randall-Sundrum ansatz [52,53]. 
Authors in [53] introduce a thickness to the brane by 
smoothing out the warp factor of a thin brane world to 
investigate the stability of a thick brane. In another ap- 
proach to derive generalized Friedmann equations, the 
four-dimensional effective brane quantities arc obtained 
by integrating the corresponding five-dimensional ones 
along the extra-dimension over the brane thickness [54] . 
These cosmological equations describing a brane of finite 
thickness interpolate between the case of an infinitely 
thick brane corresponding to the familiar Kaluza-Klein 
picture and the opposite limit of an infinitely thin brane 
giving the unconventional Friedmann equation, where the 
energy density enters quadratically. The latter case is 



then made compatible with the conventional cosmology 
at late times by introducing and fine tuning a negative 
cosmological constant in the bulk and an intrinsic pos- 
itive tension in the brane [55]. A completely different 
approach based on the gluing of a thick wall considered 
as a regular manifold to two different manifolds on both 
sides of it was first suggested in [56]. The idea behind 
this suggestion is to understand the dynamics of a local- 
ized matter distribution of any kind confined within two 
metrically different spacetimes or matter phases. Such 
a matching of three different manifolds is envisaged to 
have many diverse applications in astrophysics, early uni- 
verse, and string cosmology. Authors in [57] studied such 
a thick brane to obtain the cosmological evolution of the 
brane. Navarro and Santiago [58] considered a thick codi- 
mension 1 brane including a matter pressure component 
along the extra dimension in the energy-momentum of 
the brane. By integrating the 5D Einstein equations 
along the fifth dimension, while neglecting the parallel 
derivatives of the metric in comparison with the trans- 
verse ones, they write the equations relating the values 
of the first derivatives of the metric at the brane bound- 
ary with the integrated components of the brane energy- 
momentum tensor. These, so called matching conditions 
are then used to obtain the cosmological evolution of the 
brane which is of a non-standard type, leading to an ac- 
celerating universe for special values of the model param- 
eters. They show that when one drops the infinitesimally 
thin idealization in the modelling of the brane, gets non- 
standard cosmology on the brane. 

In Section II we make a review over the model first 
introduced by Navarro and Santiago, the cosmology of 
a thick codimension 1 brane model, its free parameters 
and modified Friedman equation which governs the back- 
ground dynamics of the universe. Then we show how 
this model can produce an accelerating universe. In Sec- 
tion III we put some constraints on the parameters of 
model by the background evolution, such as new Gold 
sample and Legacy Survey of Supernova Type la data 
[3], the combination with the position of the observed 
acoustic angular scale on CMB and the baryonic oscil- 
lation length scale. In Section IV we study the age of 
universe in the thick brane model, we also compare the 
age of the universe in this model with the age of old cos- 
mological structures in this section. Section V contains 
conclusions and discussions of this work. 

II. COSMOLOGY ON THE THICK BRANE 

In this section we make a brief review over the model 
first introduced by Navarro and Santiago to study a 
thick brane. They consider a thick brane in a Randall- 
Sundrum model with Zi symmetry. To get the cosmo- 



logical behavior of the Braneworld, the 5D metric is con- 
sidered to be as follows: 



ds 2 



n 2 (r,t)dt 2 - a 2 (r,i)dx 2 - e 2 ^dr 2 



(1) 



where r is the extra dimension and the thick brane exists 
in \r\ < e, so e is somehow the thickness of the brane. 
By the Z2 symmetry we know that the core of the brane 
is placed at r — 0. Now we should write down the 5D 
Einstein equations in the bulk and using the junction con- 
ditions derive the induced dynamics on the brane. But 
first we need to know the energy momentum tensor of 
the brane and the bulk. We assume that in the bulk 
(\r\ > e) the EMT is just that of a cosmological con- 
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And for inside the brane, a thick 



brane means we have let matter goes through the extra 
dimension. So, matter is distributed within the thickness 
of the brane and the 5D energy momentum tensor can 
be written as follows: 



1 N (brane) 



= diag(T°,T*,TyT z \T r r ) 



(2) 



where considering the matter on the brane to be a perfect 
fluid, Tq is the energy density of matter, T^,T^,T Z are 
the pressure of the matter along the normal coordinates 
and TjT is the pressure along the extra dimension. 
But, there is a problem to study this thick brane: what 
is the suitable junction condition? In the thin wall ap- 
proximation we consider the brane to be a singular hy- 
persurface with 5-function singularity in their energy- 
momentum and Einstein tensors. From the Sen-Lanczos- 
Israel (SLI) formalism we know that the difference of 
the 4D induced metric's derivatives, on both sides of the 
thin brane is given by the energy momentum tensor of 
the matter on the singular hypcrsurface. But in a thick 
brane model you can not follow the same procedure, as 
the matter is distributed along the extra dimension. So, 
what Navarro and Santiago have done is to integrate over 
Einstein equations in the \r\ < e region to obtain the suit- 
able matching conditions (they have neglected the first 
derivative of the metric along the parallel brane coordi- 
nates): 
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(Pr - P) 



(3) 
(4) 



where prime denotes differentiation with respect to r, M 3 
is the 5D fundamental mass and p, p and p r respectively 
are the 4D energy density, longitudinal (along the normal 
coordinates) and transverse 4D pressure (along the extra 
dimension), derived by integrating over the EMT along 
the thickness of the brane: 



Pr 



T*na 3 e<t>dr 



Tln^e^dr 



(6) 
(7) 



Splitting the brane EMT to the constant (A) and a time 
dependant part (with an arbitrary but assumed constant 
equation of motion), we have: 

p = A + p rn p = -A + Wp m p r = -A r + W r Pm (8) 

Now having the integrated EMT and the junction con- 
ditions we can solve the 5D Einstein equations to ob- 
tain the cosmological dynamics of the thick braneworld. 
Evaluating the rr and tr components of the 5D Einstein 
equations and keeping the terms up to the first order of 
Pm we get: 



, a a k , 
3 1- + — + —) = 

a a z a" 
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12MS 



2(A + A,.) 2 + (A + A r )(l - 3w - Aw r )p n 



Pm + 3- (1 + w) p n 

a 



WrPn 







A 

Ml 
(9) 

(10) 



where k denotes the curvature and can be k = (0, 1, —1) 
which respectively gives the spatially flat, closed and 
open universes, w and w r are the state parameters re- 
lated to the longitudinal and the transverse pressures. It 
can be easily checked that if w r = we obtain the stan- 
dard cosmology with a cosmological constant given by: 

A ?ff = 6Mf( X + X rf + IT?) U P t0 But f ° r a 

nonzero value of w r the cosmological evolution on the 
4D brane will be completely different. As it can be seen 
in the equation (10), existing a non zero pressure along 
the extra dimension violates the EMT conservation in 
the 4 dimensions: 



Pm + 3- (1 + W) Pm = WrPm + 1 

a 

One can easily integrate equation (11) to get: 



-3, 



Pm = PmOa 



(11) 



(12) 



where p m o is the present matter density. Inserting this 
result in equation (9) and integrating the equation, as- 
suming A e ff = 0, they obtain their generalized Friedman 
equation as follows: 



9 n _ 4 A + A r 2 k 

H ~ =Ca ~ + iw (1 ^ r)/,moffl "--^ 



(13) 



1 



T n °na 3 e 



b dr 



(5) 



where H = - and C is an integration constant coming 
from the extra dimension assumption. This equation 



shows that having a pressure component of EMT along 
the extra dimension is similar to having a matter with 
an effective parameter of state defined as follows: 

w + w r 



Also it can be checked from equation (9) that even 
when A e ff — we can have an accelerating universe 
{w e ff < —1/3), which means that in the case of hav- 
ing matter as the only source of the Einstein equations 
(w = 0), the evolution of the universe is an accelerat- 
ing one. So the only required condition is w r < —1/2. 
Navarro and Santiago discuss that a KK mechanism can 
produces such a negative pressure by the KK modes. So 
to have an accelerating universe we require the matter 
to be confined in the thickness of the brane which causes 
a negative w r . They have not studied the dynamics of 
the universe at early times, whether it is comparable with 
standard model of cosmology or not. As it can be checked 
easily, having w r < —1/2, all the time during the evo- 
lution of the universe, in the radiation dominated era 
(w = 1/3), equation (13) shows that the second term in 
the right hand side behaves as a b where b > —8/3, so in 
the early universe where a is very small, the first term is 
the dominated term, which corresponds to the radiation 
dominated ear a -4 . So this model is compatible with the 
standard model of cosmology at the early universe. But 
the main difference is at the late time, where the second 
term in the equation (13) is big enough and the effects 
of the brane thickness becomes important. 

In order to compare the predictions of this model with 
the observational tests, we rewrite equation (13) as a 
function of dimcnsionlcss parameters: 

H 2 (z;Q m ,n c ,w,w r ) = H^[n c (l + zf 

+fi m (l - uv)(l + zf^k - (fl tot - 1)(1 + zf] (14) 

where M 2 = = z is the rcdshift, fi m is 

the present matter density (for simplicity we discard its 
zero indice) and fl^ = O m (l — w r ) + Clc = 1 + f2fe. 

In the case of having just energy density of matter 
as the only source of the background dynamics of the 
universe, we have w = 0. For an object with the rcdshift 
of z, using the null geodesies in the FRW metric induced 
on the brane, comoving distance is obtained as: 

-no V |"fc| 

(15) 

where 

J-(x) = (x,sinx,sinhx) for k = (0, 1, —1) (16) 
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FIG. 1. Acceleration parameter (q = a/aH'^,) in the thick 
brane model as function of redshift for various values of w r . 
Here Q k = 0.0, w = 0.0 and fi ro = 0.33 

and H(z; f2 m , f^c, w r , w) is given by equation (14). To 
see how this models gives an accelerating model we study 
the behavior of the acceleration parameter (q = a/H^a). 
According to the equations (9) and (13) in this model 
accelerating parameter is derived as follows: 

q = -n c a~ 4 - -(1 + 2w r + 3w)n m a- 3( - 1+w °fri (17) 

As it can be seen in figure (1), increasing w r causes the 
universe to accelerate earlier. 

Now an interesting question that arises is: "can this 
model predict dynamics of universe?" or in another word, 
"what values of the model parameter to be consistent 
with observational tests?" 

In the forthcoming sections we will see what con- 
straints to the model described above are set by recent 
observations. 

III. OBSERVATIONAL CONSTRAINTS ON THE 
MODEL USING BACKGROUND EVOLUTION 
OF THE UNIVERSE 

In this section, at first we examine thick codimension 1 
brane model by SNIa Gold sample and supernova Legacy 
Survey data. Then to make the model parameter inter- 
vals more confined, we will combine observational results 
of SNIa distance modules with power spectrum of cosmic 
microwave background radiation and baryon acoustic os- 
cillation measured by Sloan Digital Sky survey. Table I 
shows different priors on the model parameters used in 
the likelihood analysis. 
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FIG. 2. Distance modulus of the SNIa new Gold sample 
in terms of redshift. Solid line shows the best fit values with 
the corresponding parameters of h — 0.63, Sl 7 n = 0.16+q'q3, 
tt c = 0.20t° °t and w r = -5.80+ 4 ' 80 in la level of confidence 
with Xmin/Nd.o.f = 0.92 for thick brane model. 
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FIG. 3. Distance modulus of the SNLS supernova data in 
terms of redshift. Solid line shows the best fit values with 
the corresponding parameters of h = 0.70, Q, m = 0.39+Q26' 
Sic = 0.054+° and w r = -1.80+ ' 70 in la level of confi- 
dence with Xrnin/Nd.o.f = 0.86 for thick brane model. 



TABLE I. Priors on the parameter space, used in the like- 
lihood analysis. 



Parameter 


Prior 




fltot = (1 — w r )Q, m + Sic 


1.00 


Fixed 


O.tot = (1 — Ul r )Slm + Sic 




Free 


Sim 


0.00 - 1.00 


Top hat 


Q b fi J 


0.020 ± 0.005 


Top hat (BBN) [59] 


h 




Free [60,61] 


W r 


-100.00 - 0.00 


Top hat 


W 





Fixed 


A. Supernova Type la: Gold and SNLS Samples 



The Supernova Type la experiments provided the main 
evidence of the existence of dark energy. Since 1995 two 
teams of the High-Z Supernova Search and the Super- 
nova Cosmology Project have discovered several type la 
supernovas at the high redshifts [22,62]. Recently Riess 
et al.(2004) announced the discovery of 16 type la super- 
nova with the Hubble Space Telescope. They determined 
the luminosity distance of these supernovas and with the 
previously reported algorithms, obtained a uniform 157 
Gold sample of type la supernovas [3,4,63]. Recently a 
new data set of Gold sample with smaller systematic er- 
ror containing 156 Supernova la has been released [64]. 
In this work we use this data set as new Gold sample 
SNIa. 

More recently, the SNLS collaboration released the 
first year data of its planned five-year Supernova Legacy 



Survey [65]. An important aspect to be emphasized on 
the SNLS data is that they seem to be in a better agree- 
ment with WMAP results than the Gold sample [66] . 

We compare the predictions of the thick brane model 
for apparent magnitude with new SNIa Gold sample and 
SNLS data set. The observations of supernova measure 
essentially the apparent magnitude m including redden- 
ing, K correction, etc, which are related to the (dimcn- 
sionlcss) luminosity distance, Dl-, of an object at redshift 
z through: 

m = M + 5logD L (z;fl m , fi c , w r ,w) (18) 

where 



D L (z;a 
Also 




where M is the absolute magnitude. The distance mod- 
ulus, /x, is defined as: 

fi = m — M = 5 logD^^z; f2 m , fie, w r , w) 

or 

/j, = 5logD L (z;CL m ,nc,w r ,w)+M (22) 



In order to compare the theoretical results with the 
observational data, we must compute the distance mod- 
ulus, as given by equation (21). For this purpose, the 
first step is to compute the quality of the fitting through 
the least squared fitting quantity \ 2 defined by: 



E 



\pobs(Zi) ~ fHhiZi] Hmj Q>C,W r , W, M)f 



(23) 



where cr, is the observational uncertainty in the distance 
modulus. To constrain the parameters of model, we use 
the Likelihood statistical analysis: 

C(M, Q m , Q C ,W r ,w) = M e -X 2 (M,n m ,Qc,Wr,w)/2 ( 24 ) 

where Af is a normalization factor. The parameter M is 
a nuisance parameter and should be marginalized (inte- 
grated out) leading to a new x 2 defined as: 



flk = 0.0. In this situation, the best fit values for 
the parameters of the model are Q m = O.Ol^o'Ji an d 
w r = -10.00+ 4 ' 69 with XmiJNd.o.f = 0.91 at la level 
of confidence. The corresponding value for the Hubble 
parameter at the minimized \ 2 is h = 0.63 and since 
we have already marginalized over this parameter we 
do not assign an error bar for it. The best fit values 
for the parameters of model by using SNLS supernova 
data are fi m = 0.16±^|, and w r = -4.70+ 343 with 
Xmin/Nd.o.f = 0.86 at la level of confidence. The corre- 
sponding Hubble parameter is h — 0.70. 

If we take account Qk as another free parameter, we 
would have got Vic as hoc parameter. The best fit values 
for the parameters of model by using supernova data are 



0.20: 



0.19 



and w r 



-5.80+ 4J 



0.03' 

with Xmin/Nd.o.f = 0.92 at la level of confidence. These 
values imply that ilk 
ues for the parameters of model by using SNLS super- 
nova data are O m = 0.39^;^, n c = 0.054±°°^ and 



-0.28^;^. The best fit val- 



X 2 = -21n/ C(M,Cl m ,Cl c ,w r ,w)dM (25) 

J — oo 

Using equations (23) and (25), we find: 

X 2 (0. m ,Q, c ,w r ,w) = x 2 (M = 0,n m ,Cl c ,w r ,w) 
B(n m ,rtc,w r ,w) 

c 

where 



0.053 

w r = -1.80+ - 70 with X 2 mm /N d . . f = 0.87 at la level of 
confidence. These values imply that Cl^ = —0.14 +134 . 

Figures (2) and (3) show the comparison of the theo- 
retical prediction of distance modulus by using the best 
fit values of model parameters and observational values 
from new Gold sample and SNLS supernova, respectively. 
Figures (4), (5), (6) and (7) show relative likelihood for 
ln(C/27r) (26) free parameters of thick brane model. 



B(Q m ,Cl c , w r ,w) = 

[fJ-obs(zj) - fith(zi;Vl m ,ftc,w r ,w,M = 0)] 



and 



C 



(27) 



(28) 



Equivalent to marginalization is the minimization with 
respect to M. One can show that x 2 can be expanded in 
M as [67]: 

XsNIaC^mj^C, W r ,w) = X 2 (M = 0, f2 TO , tt Cl W r , w) 

-2MB + M 2 C (29) 

which has a minimum for M = B/C: 

Xsmai^m, £Lc,w r ,w) = x 2 (M = 0,Vt m ,n c ,w r ,w) 

B(ri m ,ftc,w r ,w) 2 



c 



(30) 



Using equation (30) we can find the best fit val- 
ues of model parameters as the values that minimize 
XsNia(^»™' ^c, w r , w). In the following analysis we sup- 
pose matter domination era, w — unless stated oth- 
erwise. As a simplest case we suppose a flat universe, 



B. Combined analysis: SNIa+CMB+SDSS 

Before last scattering, the photons and baryons are 
tightly coupled by Compton scattering and behave as 
a fluid. The oscillations of this fluid, occurring as a 
result of the balance between the gravitational interac- 
tions and the photon pressure, lead to the familiar spec- 
trum of peaks and troughs in the averaged temperature 
anisotropy spectrum which we measure today. The odd 
peaks correspond to maximum compression of the fluid, 
the even ones to rarefaction [68] . In an idealized model of 
the fluid, there is an analytic relation for the location of 
the ra-th peak: l m ps ml a [69,70] where I a is the acoustic 
scale which may be calculated analytically and depends 
on both pre- and post-recombination physics as well as 
the geometry of the universe. The acoustic scale corre- 
sponds to the Jeans length of photon-baryon structures 
at the last scattering surface some ~ 379 Kyr after the 
Big Bang [13] . The apparent angular size of acoustic peak 
can be obtained by dividing the comoving size of sound 
horizon at the decoupling epoch r s (zdec) by the comoving 
distance of observer to the last scattering surface r(zdec)'- 



Tj_ = r s {z dec ) 
I A r(z dec ) 



(31) 



The size of sound horizon at the numerator of equation 
(31) corresponds to the distance that a perturbation of 



pressure can travel from the beginning of universe up to 
the last scattering surface and is given by: 



v s (z')dz' 
H(z')/H 



(32) 



where v s (z)~ 2 = 3 + 9/4 x Pb{z) / p ra d(z) is the sound 
velocity in the unit of speed of light from the big bang 
up to the last scattering surface [24,69] and the rcdshift 
of the last scattering surface, z<; ec , is given by [69]: 



z dec = 1048 [1 + 0.00124(^)-°- 738 ] [1 + .91 (^m) 92 
,9i = 0.0783(cj 6 )-°- 238 Fl 1 5nt: '"^ ' 763 ' 
.92 = 0.560 [l + 21.1(w 6 ) 1 ' 81 ] 



[1 + 39.5(wh) c 

, -l 



(33) 



where uj m = D, m h 2 , p ra d is the radiation density and tOb = 
Qbh 2 (fib is the present baryonic density). Changing the 
parameters of the model can change the size of apparent 
acoustic peak and subsequently the position of I a = k/9a 
in the power spectrum of temperature fluctuations on 
CMB. The simple relation l m « ml a however does not 
hold very well for the first peak although it is better 
for higher peaks. Driving effects from the decay of the 
gravitational potential as well as contributions from the 
Doppler shift of the oscillating fluid introduce a shift in 
the spectrum. A good paramcterizations for the location 
of the peaks and troughs is given by [25,70] 



l m = U(m - <j) m ) 



(34) 



where <p m is phase shift determined predominantly by 
pre-recombination physics, and are independent of the 
geometry of the Universe. The location of first acoustic 
peak can be determined in model by equation (34) with 
</>i(w m , o->b) ~ 0.27 [25,70]. Instead of the peak locations 
of power spectrum of CMB, one can use another model 
independent parameter which is so-called shift parameter 
TZ as: 



TZ oc 



flat 



h 



(35) 



where l[ lat corresponds to flat pure-CDM model with 
f2 m = 1.0 and the same uj m and cob as the original model. 
It is easily shown that shift parameter is as follows [71] : 



TZ = 



- Dh{zdec\ O m ,Oc,W r ,H)) 
(1 + Z dec ) 



(36) 



The observational results of CMB experiments corre- 
spond to a shift parameter of TZ = 1.716 ± 0.062 (given 
by WMAP, CBI, ACBAR) [13,72]. One of the advan- 
tages of using the parameter TZ is that it is independent 
of Hubble constant. In order to put constraint on the 
model from CMB, we compare the observed shift param- 
eter with that of model using likelihood statistic as [71]: 



C ^ g~XcMB/ 2 



where 



Xcmb 



[7^ bs — TZtxX 

2 

°CMB 



(37) 



(38) 



where TZth and TZ \, S are the theoretical shift parameter, 
determined using equation (36), and the observed one, 
respectively. 

The large scale correlation function measured from 
46, 748 Luminous Red Galaxies (LRG) spectroscopic 
sample of the SDSS (Sloan Digital Sky Survey) includes 
a clear peak at about 100 Mpc h~ x [73]. This peak was 
identified with the expanding spherical wave of baryonic 
perturbations originating from acoustic oscillations at re- 
combination. The comoving scale of this shell at recom- 
bination is about 150Mpc in radius. In other words this 
peak has an excellent match to the predicted shape and 
the location of the imprint of the recombination-epoch 
acoustic oscillation on the low-redshift clustering matter 
[73] . A dimcnsionless and independent of Hq version of 
SDSS observational parameter is: 



A = D v (z sdss ) 



^JVl m H 2 

Zsdss 

H D 2 L (z sdss ; n m ,tt c ,w r , w) 



_H(z sdss -MmMc,w r ,w)z 2 dss (l + z sdss ) 2 



1/3 



(39) 



where Dv(z sdss ) is characteristic distance scale of the 
survey with mean redshift z s dss [73-75]. We use the 
robust constraint on the thick branc model using the 
value of A = 0.469 ± 0.017 from the LRG observation 
at z s dss = 0.35 [73]. This observation permits the addi- 
tion of one more term in the % 2 of equations (30) and (38) 
to be minimized with respect to H(z) model parameters. 
This term is: 



XsDSS 



[A Q hs ^ AthY 



(40) 



SDSS 



This is the third observational constraint for our analysis. 

In what follows we perform a combined analysis of 
SNIa, CMB and SDSS to constrain the parameters of 
the thick branc model by minimizing the combined x 2 = 
XsNia + Xcmb + Xsdss- Thc bcst values of the model 
parameters from the fitting with the corresponding er- 
ror bars from the likelihood function marginalizing over 
the Hubble parameter in the multidimensional parameter 
space in flat universe and domination matter epoch are: 



1.92±2il at la confidence 



fi m = 0.331^3 and m 



level with Xmin/Nd.o.f = 1-02. The Hubble parameter 
corresponding to the minimum value of x 2 is ft, = 0.64. 
Here we obtain an age of 15.62^3'g^ Gyr for the uni- 
verse (see section IV for more details). Using thc SNLS 
data, the best fit values of model parameters are: fl m = 



TABLE II. The best values for the parameters of a thick 
brane model with the corresponding age for the universe from 
fitting with SNIa from new Gold sample and SNLS data, 
SNIa+CMB and SNIa+CMB+SDSS experiments at one and 
two a confidence level. Here we suppose Qk = 0.0 and 



w = 0.0. 



Observation 




Wr 


Age (Gyr) 


SNIa(new Gold) 


01+ 012 
01+°' 33 


-10.00+ 469 
— 10.00 +98 ' 31 


8.06+ 314 


SNIa(new Gold)+CMB 


43+™! 

0.08 


-1 24+ 022 

■"■•^-O-SO 

-1 24+ ,- 39 

'■• en — 0.47 


15.2ll|;« 


SNIa(new Gold) + 
CMB+SDSS 


n Qo+0.02 
u.oo_ 03 

n OO+0.05 
u.oo_ 05 


-1 q 2 + - 18 

-1 Q2+ 037 


15-6211;^ 


SNIa (SNLS) 


U.1D_ 16 

16+ 035 


-4.70+ 3 - 43 
-4.70+ 373 


12.64+ 874 


SNIa(SNLS)+CMB 


n qo+0.05 
u - oo -0.05 

u -°°-0.09 


-i.9HS:S 
-i.9il!:3 


1 4 sn+ 4 - 49 

r^.ou_ 4 86 


SNIa(SNLS)+ 
CMB+SDSS 


n os +o.02 

U.^8_ .02 


-2 47+ - 27 

AAi -0.30 

o 47+ - 52 


1 c; no+4.04 
10.ZO_ 4 07 



0.28Io;o2 and w r = -2.47l2:3o at lcr confidence level with 
Xmin/Nd.o.f = 0.86. Age of universe calculating with the 
best fit parameters is 15. 231^07 ( see nex t section). Ta- 
ble II indicates the best fit values for the cosmological 
parameters with one and two a level of confidence. In 
general case, the following values for the free parameters 
maximize likelihood probability: Sl Tn = 0.311q'q2, Sic = 
0.05t£oi and w r = -1.4Q±%%% at lcr confidence level 
stats n k = +0.2ll°;°|. SNLS SNIa+CMB+SDSS give: 
Q m = 0.28l°;^, n c = 0.037lg:gg2 and «v = -2.05±°;^ 
at lcr confidence level demonstrate fl k = +0.1ll : 07 . Ta- 
bles III and IV give the best fit values for free parameters 
and age of universe computing with these values. Joint 
confidence intervals in free parameter spaces are shown 
in figures (8), (9), (10), (11), (12) and (13). 



TABLE III. The best fit values for the parameters of the 
model using SNIa from new Gold sample and SNLS data, 
SNIa+CMB and SNIa+CMB+SDSS experiments at one and 
two a confidence level. Here we suppose w = 0.0. 



Observation 




fie- 




SNIa(new Gold) 


16+ 084 

u - lo -0.03 

16+ 084 

u - lo -0.08 


ri 20+ 004 


-5.80+ 4 - 80 
-5.80+ 4 - 90 


SNIa(new Gold)+CMB 


n S3 +0 - 10 
n S3+ 016 


03+ 02 

U.UO_ 02 

o.o3i2:2 3 3 


-1 io +0 - 20 

i - lu -0.20 

-1 io +0 - 40 


SNIa(new Gold)+ 
CMB+SDSS 


S1+ 02 

u ' o± -0.02 
u - o± -0.02 


o.o5±2:2i 
o.o5l2:2i 


_i 4n+ - 20 
-l 4n+ - 20 

l - w -0.20 


SNIa (SNLS) 


SQ+ - 61 

u -° 3 -0.30 


n n54+ 056 

U.UO*_ .053 

n n54+ 086 

U.U04_o. 053 


-1.80+ - 70 

_ L80 +0.80 


SNIa(SNLS)+CMB 


42+ 09 

"•^-0.23 

42+ 12 

"•^-0.33 


n 027+ 013 

u.u^l -0.023 

n 027+ ' 018 

u.u^l -0.026 


-1 60 +a3 ° 

L - ou -1.05 

-i-6o«: 4 o2 


SNIa(SNLS)+ 
CMB+SDSS 


0.28lg;S 
0.28l ;2I 


U,UJ ' -0.004 

n o 7 +0.008 
U.UO( _ o8 


-2.0512:^ 
-2.05«: 3 ^ 



IV. AGE OF UNIVERSE 

The age of universe integrated from the big bang up 
to now in terms of free parameters of thick brane model 
is given by: 

t (Q m ,Qc,w r ,w) = I dt 
Jo 

1 ( r, r f°° dz ' H o \ / x 

'iwmfi^l ottwm) (41) 

Figure (14) shows the dependence of H t (Hubble pa- 
rameter times the age of universe) on Sl m and w r for a 
flat universe. Obviously increasing Sl m and w r result in a 
longer and shorter age for the universe, respectively. As 
a matter of fact, according to the equation (14), f2 m |izv| 
behaves as dark energy in the ACDM scenario and w r 
has the same role as w in the ACDM (see figures (15)). 



TABLE IV. The best values for the curvature of a thick 
brane model with the corresponding age for the universe from 
fitting with SNIa from new Gold sample and SNLS data, 
SNIa+CMB and SNIa+CMB+SDSS experiments at one and 
two a confidence level. Here we suppose w = 0.0. 



f 1 Kc^srT 7"Q "f" 1 /"111 


o, 




SNIa(new Gold) 


-outfit 


12 54+8.85 


SNIa(new Gold)+CMB 


+0.14±°;« 


15.3111;°? 


SNIa(new Gold) + 
CMB+SDSS 


+0.21±°;£ 


14 7 2+0'« 

' z -0.48 


SNIa (SNLS) 


-0.14+ 134 


13.50+ 295 


SNIa(SNLS)+CMB 


u - iu -0.74 




SNIa(SNLS)+ 
CMB+SDSS 


+0 11+ - 10 


14.18lg;g 



The "age crisis" is one the main reasons of the accel- 
eration phase of the universe. The problem is that the 
universe's age in the Cold Dark Matter (CDM) universe 
is less than the age of old stars in it. Studies on the old 
stars [76] suggest an age of 13+2 Gyr for the universe. 
Richer et. al. [77] and Hasen et. al. [78] also proposed 
an age of 12.7 ± 0.7 Gyr, using the white dwarf cool- 
ing sequence method (for full review of the cosmic age 
see [13]). To do another consistency test, we compare 
the age of universe derived from this model with the age 
of old stars and Old High Rcdshift Galaxies (OHRG) in 
various rcdshifts. Table IV shows that age of the uni- 
verse from the combined analysis of SNIa+CMB+SDSS 
is 14.72 j£g Gyr and 14.18±g || Gyr for new Gold sam- 
ple and SNLS data, respectively, while ACDM implies 
13.7 ± 0.2Gyr [13]. These values are in agreement with 
the age of old stars [76] . 

Here we consider three OHRG for comparison with 
the thick brane model, namely the LBDS 53W091, a 
3.5-Gyr old radio galaxy at z = 1.55 [79], the LBDS 
53W069 a 4.0-Gyr old radio galaxy at z = 1.43 [80] and 
a quasar, APM 08279 + 5255 at z = 3.91 with an age 
of t = 2.1+Q ^Gyr [81]. The latter has once again led to 
the " age crisis" . An interesting point about this quasar 
is that it cannot be accommodated in the ACDM model 
[82]. In order to quantify the age-consistency test we 
introduce the expression r as: 



TABLE V. The value of r for three high redshift objects, 
using the parameters of the model derived from fitting with 
the observations. Here we imagined flat universe and w = 0.0. 



Observation 


LBDS 
53W069 
z = 1.43 


LBDS 
53W091 
z = 1.55 


APM 
08279 + 5255 
z = 3.91 


SNIa (new Gold) 


u - od -0.02 


S6+ 013 

u -°°-0.02 


16+ 007 


SNIa(new Gold)+CMB 
+SDSS 


i 9 7 +0.54 
L -^' -0.52 


i - JD -0.57 


u - ' J -0.37 


SNIa(new Gold)+CMB 
+SDSS+LSS 


i Qo+0.48 
1 -°°-0 .53 




n 7 5+0.36 
u - ' J -0.38 


SNIa (SNLS) 


n 85+ - 64 


o.89lg:£ 


u -^ z -0.42 


SNIa(SNLS)+CMB 
+SDSS 


1 1Q+0-61 
1 - ±y -0.64 


i 9 7+0.65 
L -^' -0.69 


0.67iS;S 


SNIa(SNLS)+CMB 
+SDSS+LSS 


1 36+ 62 

l.OU_ 62 


i - 4O -0.67 


0.78l° ;« 


TABLE VI. The value of r for three high redshift objects, 
using the parameters of the model derived from fitting with 
the observations. Here we chose w — 0.0. 


Observation 


LBDS 
53W069 
z = 1.43 


LBDS 
53W091 
z = 1.55 


APM 
08279 + 5255 
z = 3.91 


SNIa (new Gold) 


n 7 c+0.44 
U - '°-0.75 


o.78ig:« 




SNIa(new Gold)+CMB 
+SDSS 


o 7 +0.23 
L -°' -0.41 




0.83i° ; 24 


SNIa(new Gold)+CMB 
+SDSS+LSS 


-1 ,0+0.07 
1 -- LJ -0.08 


1 1Q+0.07 


0.65i° ;i 4 


SNIa (SNLS) 


i m +0 - 34 

J — LU -1.10 


i - LD -1.16 


u -°°-0.58 


SNIa(SNLS)+CMB 
+SDSS 


1 or+0.25 


1 44+0-27 


0-80i°; 22 


SNIa(SNLS)+CMB 
+SDSS+LSS 


i 7+0.05 
i ' 1 '-0.05 


1 94+0.06 


67+°' 14 
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FIG. 4. Marginalized likelihood functions of two parame- 
ters of model (Q m and w r ). The solid line corresponds to the 
likelihood function of fitting the model with SNIa data (new 
Gold sample), the dashdot line with the joint SNIa+CMB 
data and dashed line corresponds to SNIa+CMB+SDSS. 
The intersections of the curves with the horizontal solid and 
dashed lines give the bounds with la and 2a level of confi- 
dence respectively. Here = 0.0 and w — 0.0. 

_ t(z;n m ,Q, c ,w r ,w) _ t(z;il m ,rtc,w r ,w)H 

"Obs vohg-LJ-Q 

where t(z) is the age of universe, obtained from the equa- 
tion (41) and t (, s is an estimation for the age of old 
cosmological object. In order to have a compatible age 
for the universe we should have r > 1. Tables V and 
VI report the value of r for three mentioned OHRG with 
various observations. We see that the parameters of thick 
brane model from the combined observations provide a 
compatible age for the universe, compared to the age 
of old objects, also in addition SNLS data result in a 
shorter age for the universe. Once again for the thick 
brane model, APM 08279 + 5255 at z = 3.91 has a longer 
age than the universe but gives better results than most 
cosmological models investigated before [14,15,82]. 




FIG. 5. Marginalized likelihood functions of two param- 
eters of thick brane model (fi m and w r ). The solid 
line corresponds to the likelihood function of fitting the 
model with SNIa data (SNLS), the dashdot line with the 
joint SNIa+CMB data and dashed line corresponds to 
SNIa+CMB+SDSS. The intersections of the curves with the 
horizontal solid and dashed lines give the bounds with la 
and 2a level of confidence respectively. Here Qt = 0.0 and 
w = 0.0. 

V. CONCLUSIONS AND DISCUSSIONS 

From observational point of view, it has been possi- 
ble to compare theoretical model with the observational 
results. 

We explored the consistency of a thick codimcnsion 1 
brane model with the implication of up-to-date luminos- 
ity of supernova type la observed by two independent 
groups, new Gold sample and SNLS data set, acous- 
tic peak in the cosmic microwave background anisotropy 
power spectrum and baryon acoustic oscillation measured 
by Sloan Digital Sky Survey. 

In this scenario, universe is supposed to be dominated 
by pressureless cold dark matter, which penetrates to the 
extra dimension, this leads to an acceleration epoch for 
the universe. 
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FIG. 6. Marginalized likelihood functions of three pa- 
rameters of model (fi m , w r and fic). The solid line cor- 
responds to the likelihood function of fitting the model 
with SNIa data (new Gold sample), the dashdot line with 
the joint SNIa+CMB data and dashed line corresponds to 
SNIa+CMB+SDSS. The intersections of the curves with the 
horizontal solid and dashed lines give the bounds with la and 
2a level of confidence respectively. Here we take w = 0.0. 




FIG. 7. Marginalized likelihood functions of three param- 
eters of model (Sl m , w r and Sic). The solid line corresponds 
to the likelihood function of fitting the model with SNIa data 
(SNLS), the dashdot line with the joint SNIa+CMB data and 
dashed line corresponds to SNIa+CMB+SDSS. The intersec- 
tions of the curves with the horizontal solid and dashed lines 
give the bounds with la and 2a level of confidence respec- 
tively. Here we take w = 0.0. 




FIG. 8. Joint confidence intervals of S7 m and w r , fitted with 
SNIa new Gold sample+CMB+SDSS. Solid line, dashed line 
and long dashed line correspond to 3a, 2a and la level of 
confidence, respectively. Here Qt = 0.0 and w = 0.0. 
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FIG. 9. Joint confidence intervals of fl m and w r , fitted with 
SNIa SNLS+CMB+SDSS. Solid line, dashed line and long 
dashed line correspond to 3<r, 2a and la level of confidence, 
respectively. Here fi^ = 0.0 and w — 0.0. 



In this work we have been interested in matter dom- 
inated era for the universe. So we imagined w = 0.0 
as a prior through this paper. The best parameters 
obtained from the fitting with the new Gold sample 
data combined with CMB and SDSS observations are: 
n rn = 0.3lt°;°i n c = 0.05±°;°* and w r = -1.40±°|° 
at la confidence level states spatially open universe 
with n k = +0.21±g;gf. SNLS SNIa+CMB+SDSS give: 
fl m = 0.28±£^, Q c = 0.037±g ggi| and w r = -2.05t 015 



0.15 



O.ll+o 1 ;^. The 



at la confidence level demonstrate ilk 
well-known ACDM model implying -0.06 < tt k < +0.02 
[13] and some other interesting models such as Dvali- 
Gabadadze-Porrati (DGP) which states fl k = O.Ol^og 
and Q k = 0.011^4 

using Gold sample and SNLS data, 
respectively [83,84], show a contradiction with our re- 
sults. In fact having a spatially open universe is ruled 
out in many models comparing to the observations while 
in this thick brane model we find out it is not possible to 
have a spatially flat universe, according to the recent ob- 
servational tests. The value of w r given by observational 
constraints is negative. This shows that instead of dark 
energy to accelerate the universe, we have a strange ef- 
fect of matter through the extra dimension with negative 
pressure. 

We also performed the age test, comparing the age 
of old stars and old high redshift galaxies with the age 
derived from this model. From the best fit parameters 
of the model using new Gold sample and SNLS SNIa, 
respectively, we obtained an age of 14.72^Q 4g Gyr and 
14. 18^0^9 Gyr, for the universe. These results are in 
agreement with the age of the old stars. The age of uni- 
verse in this model is larger than what is given in the 
other models [13-15,84] 

To check the age crisis in this model we chose two high 
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FIG. 10. Joint confidence intervals of Q m and w r , fitted 
with SNIa new Gold sample+CMB+SDSS. Solid line, dashed 
line and long dashed line correspond to 3a, 2a and la level 
of confidence, respectively. Here we supposed w = 0.0. 



redshift radio galaxies at z = 1.55 and z — 1.43 with 
a quasar at z = 3.91. Two first objects were consistent 
with the age of universe, i.e., they were younger than 
the universe while the third one was not but gave better 
result than ACDM and a class of Quintessence model 
[14,15]. 
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